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High temperature QCD with three flavors of improved staggered quarks * 
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We present an update of our study of high temperature QCD with three flavors of quarks, using a Symanzik 
improved gauge action and the Asqtad staggered quark action. Simulations are being carried out on lattices with 
Nt = 4, 6 and 8 for the case of three degenerate quarks with masses less than or equal to the strange quark 
mass, m 3 , and on lattices with Nt — 6 and 8 for degenerate up and down quarks with masses in the range 
0.2m s < m Ut d < 0.6m s , and the strange quark fixed near its physical value. We also report on first computations 
of quark number susceptibilities with the Asqtad action. These susceptibilities are of interest because they can be 
related to event-by-event fluctuations in heavy ion collision experiments. Use of the improved quark action leads 
to a substantial reduction in lattice artifacts. This can be seen already for free fermions and carries over into our 
results for QCD. 



We are studying high temperature QCD with 
three flavors of improved staggered quarks 
using a one-loop Symanzik improved gauge ac- 
tion and the Asqtad quark action [^|j3|. Both 
the gauge and quark actions have all lattice arti- 
facts removed through order a 2 (a is the lattice 
spacing) at tree level, and are tadpole improved. 
Thus, the leading order finite lattice spacing arti- 



facts are of order a. 



We consider two cases: 



1) all three quarks have the same mass m q ; and 

2) the two lightest mass quarks have equal mass 
m Uj d and the mass of the third quark is fixed at 
that of the strange quark m s . We refer to these 
cases as Nf — 3 and Nf = 2 + 1, respectively. 

Because it contains the Naik term the Asq- 
tad action has a significantly better dispersion re- 
lation for the quarks than the standard Kogut- 
Susskind and Wilson quark actions, which helps 
decrease lattice artifacts in the energy and pres- 
sure [0 , and, as we shall see, for the quark number 
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susceptibilities. The Asqtad action also exhibits 
excellent scaling properties ||, and significantly 
better flavor symmetry properties than the con- 
ventional Kogut-Susskind action We esti- 
mate that for lattices with eight to ten time slices, 
the kaon will be heavier than the heaviest non- 
Goldstone pion in the neighborhood of the finite 
temperature transition or crossover. 

We have attempted to vary the temperature 
while keeping all other physical quantities con- 
stant. To this end we have performed a set 
of spectrum calculations at a — 0.13 fm and 
0.20 fm for three equal mass quarks. Here, and 
throughout this work, we determine the lattice 
spacing from the static QQ potential, and ex- 
press dimensionful quantities in terms of r\ de- 
fined by r\FQq{r\) = 1, which leads to r\ w 
0.34 fm. We determine m s from the requirement 
that m rhs jra^ w 0.673. We have performed spec- 
trum calculations at both lattice spacings with 
m„ = m s , 0.6 m s and 0.4 m s , with an additional 



2 



calculation at 0.2 m s in progress. We would like 
to carry out our thermodynamics studies with 
three equal mass quarks for m^ as /m<^ fixed, but 
this quantity will, of course, vary slightly with 
lattice spacing if we keep m q /m s fixed. So, at 
a = 0.2 fm we perform linear interpolations of 
m^ ss and in the quark mass to determine the 
precise values of m q for which m r)gg /m^ takes on 
the values found at a = 0.13 fm. Then to de- 
termine the values of am q and a for thermody- 
namics studies with 0.13 fm < a < 0.20 fm, we 
perform linear interpolations of the logarithms of 
these quantities in the gauge coupling 10/g 2 . 

Our approach for thermodynamics studies with 
mu,d < Tn s is quite similar. In this case we wish 
to vary the temperature keeping both m v /m p 
and ro, si /m^, fixed. We have therefore performed 
a set of spectrum calculations at a — 0.13 fm 
and 0.20 fm with the mass of the light quarks, 
fn,u,d — 0.6 77i s , 0.4m s and 0.2 m s . An addi- 
tional spectrum calculation with m Ui d = 0.1 m s 
is in progress. In these calculations the mass of 
the heavy quark is fixed at m s , as determined 
from spectrum calculations with three equal mass 
quarks. In our spectrum runs at a — 0.13 fm |(| 
we found that varied by less than 2% and 

rrirf, by less than 1% for 0.2m s < m u ,d < Ti s and 
the heavy quark mass held fixed at m s . So, the 
neglect of the dependence of m, M and on m u ,d 
is well justified. In these studies, we performed 
linear interpolations of m„ and m p at a = 0.20 fm 
to determine the values of m Ui d for which m n /m p 
takes on the values found at a = 0.13 fm. Then 
for 0.13 fm < a < 0.20 fm we perform linear 
interpolations of the logarithms of a, am s , and 
arn Ut d in the gauge coupling along lines of con- 
stant m^/nip and rn Vss /m ( j,. 

In all figures below we give the values oim q /m s 
for Nf = 3 and m Ui d/m s for Nf — 2 + 1 at a = 
0.13 fm. The corresponding values of m^Jm^ 
and m^/rrip are given in Table |. 

For three equal mass quarks, Nf = 3, we have 
carried out thermodynamics studies on lattices 
with four, six and eight times slices, and aspect 
ratio N a /Nt = 2. Here N s and N t are the spatial 
and temporal dimensions of the lattice in units 
of the lattice spacing. The spectrum calculations 
and interpolations described above allowed us to 



Table 1 

In the first column we show the value of m q /m s 
at lattice spacing 0.13 fm, which produced the 
m riaa /m ( j > ratio shown in the second column for 
spectrum calculations with three equal mass 
quarks. In the third column we give the value of 
mu^d/rris which produced the m w /m p ratio shown 
in the fourth column for spectrum calculations 
with two equal mass light quarks, and the mass 
of the heavy quark fixed at m s . 



N f = 3 


N f = 


2 + 1 


m q /m s 




m Ut d/m s 


m w /m p 


1.0 


0.673 


1.0 


0.673 


0.6 


0.583 


0.6 


0.582 


0.4 


0.504 


0.4 


0.509 






0.2 


0.392 



determine the values of the quark mass am q that 
keep m^ 3s /m0 fixed as the gauge coupling is var- 
ied. They also enabled us to determine the value 
of the lattice spacing and, therefore, the temper- 
ature for each run. 

In Fig. [l] we plot the real part of the Polyakov 
loop as a function of temperature on 16 3 x 8 lat- 
tices. The Polyakov loop shows a crossover from 
confined behavior at low temperature to decon- 
fined behavior at high temperature, but the tran- 
sition is not particularly sharp. The insensitiv- 
ity of the Polyakov loop to the quark mass at 
fixed temperature or lattice spacing is perhaps 
not surprising. We have determined the lattice 
spacing from the heavy quark potential, and in 
our spectrum runs we adjusted the coupling con- 
stant to keep the lattice spacing fixed as the quark 
mass is varied. Since the Polyakov loop, like the 
heavy quark potential, is determined from mea- 
surements of purely gluonic operators, our proce- 
dure is likely to minimize the dependence of the 
Polyakov loop on the quark mass. 

In Fig. |^ we show the chiral order parameter, 
iptp as a function of temperature for m q — m s , 
0.6 m s and 0.4 m s on 16 3 x 8 lattices. The black 
bursts in this figure are linear extrapolations of 
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Figure 1. Real part of the Polyakov loop on 16 3 x 
8 lattices for three degenerate flavors of quarks. 
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Figure 2. The chiral order parameter, ifnp, 16 3 x 8 
lattices. The black bursts are linear extrapola- 
tions in the quark mass to m q — for fixed tem- 
perature. 
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Figure 3. Real part of the Polyakov loop on 16 3 x 
8 lattices for two light and one heavy quark. 



tp%/} in the quark mass to m q — for fixed tem- 
perature. All of these results are indicative of 
a crossover at the quark masses studied so far, 
rather than a bona fide phase transition. This re- 
sult is consistent with earlier studies which found 
phase transitions only for smaller values of the 
quark mass than we have studied to date 0. 
Fig. H suggests that there is unlikely to be a phase 
transition for temperatures above 180 MeV, but 
one could occur at or below that value. 

The Nf = 2 + 1 thermodynamics studies were 
carried out on 12 3 x 6 and 16 3 x 8 lattices. In 
this phase of our work we performed simulations 
with two degenerate light quarks and the heavy 
quark mass held equal to that of the strange 
quark. The spectrum calculations and interpo- 
lations described above allowed us to determine 
the values of the heavy and light quark masses 
needed to keep m^/mp and rn r]B J'm§ fixed as the 
gauge coupling was varied. They also allowed us 
to determine the lattice spacing, and therefore the 
temperature for each run. 

In Fig. H we plot the real part of the Polyakov 
loop as a function of temperature on 16 3 x 8 lat- 
tices for the three values of m Uj( j studied to date. 
As in the Nf = 3, study we observe a crossover 
from confined to deconfined behavior, rather than 
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Figure 4. The chiral order parameter, ipip, on 
16 3 x 8 lattices The black bursts are linear ex- 
trapolations in the quark mass to m q — for 
fixed temperature. 
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Figure 5. The quark number susceptibility for 
free Asqtad and standard staggered fermions as 
function of temporal lattice size Nt ■ 



a sharp transition, and little dependence on the 
light quark mass. 

In Fig. H we plot the chiral order parameter as 
a function of temperature on 16 3 x 8 lattices. The 
green, blue and red plotting symbols are data for 
Tn Ut d — 0.6 m s , 0Am s and 0.2 m s respectively, 
and the black bursts are linear extrapolations in 
vnu,d for fixed temperatures. As in the case of 
Nf — 3, all of these results are indicative of 
a crossover at the quark masses studied so far, 
rather than a bona fide phase transition. Fig. ^ 
suggests that for the heavy quark mass equal to 
that of the strange quark, there is unlikely to be 
a phase transtion for T > 170 MeV. Simulations 
at smaller values of m u ,d will be required to de- 
termine the precise location and the nature of the 
expected finite temperature transition. We plan 
to carry out such simulations during the coming 
year. 

We have also measured quark number suscep- 
tibilities They can be related to event-by- 
event fluctuations in heavy ion collisions fllPf by 
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Figure 6. The triplet quark number susceptibility 
for Nf = 3 on 12 3 x 6 lattices. The black line on 
the right of the figure indicates the free quark 
value in the continuum, and the red line the free 
quark value on the finite lattices used. 
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Figure 7. The triplet quark number susceptibility 
for N f = 2 + 1 on 16 3 x 8 lattices. The black line 
on the right of the figure indicates the free quark 
value in the continuum, and the red line the free 
quark value on the finite lattices used. 



the fluctuation-dissipation theorem 
T d 2 \ogZ 



(5Q 2 ) 



V s dfi 2 Q 



Xq{T,hq = 0) , 



(1) 



with fig the chemical potential for the conserved 
charge Q, e.g. strangeness. We introduce a chem- 
ical potential fi into the Asqtad action by 

p x , y (tx) = Jp svatlal + (2) 
m (x) [u<P(x)erS x+6ty - C/f )f (x)e^S xy+() 

+U^ L \x)e^S x+ ^ y - U^Hx)e~^8 x ^_ , 

where Uq is the time-like "fat-link" and Uq L ^ 
the time-like "long-link" including their tadpole 
improved coefficients. 
Defining 



Xij 



T d 2 logZ 



V s dnidfXj 



(3) 



we computed the light- 2-flavor singlet 

Xsing 

= 2 X uu + 2 X ud (4) 
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Figure 8. The triplet quark number susceptibility 
for Nf = 3 with quarks of mass m q = 0.4m s , on 
8 3 x 4, 12 3 x 6 and 16 3 x 8 lattices. The black 
line on the right of the figure indicates the value 
for free quarks in the continuum, and the colored 
lines the value for free quarks on the finite lattices 
used. 



and triplet quark number susceptibility 

Xtrip — ^Xuu ^Xud ? 



(5) 



using stochastic estimates with 100 random Gaus- 
sian sources. For Nf = 2 + 1 we also com- 
puted the strange quark number susceptibility 

Xstrange — Xss • 

For free fermions, the Asqtad action, because 
of the Naik term, reduces lattice artifacts signifi- 
cantly compared with the standard staggered ac- 
tion as demonstrated in Fig. ^. 

In Figs. |^ and we show the triplet suscepti- 
bility for the three quark masses we have studied 
in our Nf = 3 simulations on 12 3 x 6 lattices, 
and for the three values of m Ut d we have studied 
in our Nf = 2 + 1 simulations on 16 3 x 8 lat- 
tices, respectively. In Figs. || and [| we show the 
triplet susceptibility for fixed quark mass for the 
different values of Nt used. In Fig. ^ we also show 
2 x Xss/T 2 , a quantity which should become equal 
to Xtrip /T 2 in the very high temperature limit. 
We observe that the quark number susceptibil- 
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Figure 9. The triplet and strange quark number 
susceptibilities for Nj = 2 + 1 with m U: d — 0.2 m s 
on 12 3 x 6 and 16 3 x 8 lattices. The black line on 
the right of the figure indicates the value for free 
quarks in the continuum, and the colored lines the 
value for free quarks on the finite lattices used. 

ity provides a clean signal for the crossover from 
the hadronic to the quark-gluon plasma phase. 
The close agreement between the N t = 6 and 8 
results in these figures is an indication of the ex- 
cellent scaling properties in lattice spacing of the 
Asqtad action. 

The singlet susceptibility looks very similar to 
the triplet susceptibility shown in Figs. || and ^. 
We show in Fig.|l(] the difference between singlet 
and triplet susceptibility for Nf — 3 at the light- 
est quark mass considered so far, m q = 0.4 m s . 
The two are clearly different in the crossover re- 
gion. 

This work is supported by the US National Sci- 
ence Foundation and Department of Energy and 
used computer resources at Florida State Univer- 
sity (SP), NCSA, NERSC, NPACI, FNAL, and 
the University of Utah (CHPC). 

REFERENCES 

1. Earlier reports on this work can be found in 
C. Bernard et al. (The MILC collaboration), 



o.o 



E ' 



X 



-0.2 — 




100 



200 300 

T (MeV) 



400 



Figure 10. The difference between singlet and 
triplet quark number susceptibility for Nf = 3 
with quarks of mass m q — OA m s , on 8 3 x4, 12 3 x6 
and 16 3 x 8 lattices. 

Nucl. Phys. A 702 (2002) 140; Nucl. Phys. B 
(Proc. Suppl) 106 (2002) 429; 

2. K. Orginos, D. Toussaint and R.L. Sugar, 
Phys. Rev. D 60 (1999) 054503; Nucl. Phys. 
(Proc. Suppl.) 83 (2000) 878. 

3. CP. Lepage, Phys. Rev. D 59 (1999) 074501. 

4. S. Naik, Nucl. Phys. B316 (1989) 238. 

5. C. Bernard et al. (The MILC collaboration), 
Phys. Rev. D 61 (2000) 111502. 

6. C. Bernard et al. (The MILC collaboration), 
Phys. Rev. D 64 (2001) 054506. 

7. S. Aoki et al. (JLQCD), Nucl. Phys. B (Proc. 
Suppl) 73 (1999) 459; F. Karsch, E. Laer- 
mann and Ch. Schmidt, Phys. Lett. B520 
(2001) 41. 

8. S. Gottlieb et al. Phys. Rev. Lett. 59 (1987) 
2247. 

9. R.V. Gavai and S. Gupta, Phys. Rev. D 
64 (2001) 074506; Phys. Rev. D 65, (2002) 
094515; R.V. Gavai, S. Gupta and P. Majum- 
dar, Phys. Rev. D 65 (2002) 054506. 

10. B. Miiller, Nucl. Phys. A 702 (2002) 281; V. 
Koch, M. Bleicher and S. Jeon, Nucl. Phys. 
A 702 (2002) 291; R.V. Gavai, Nucl. Phys. A 
702 (2002) 299. 



